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Abstract 

The recently invented cylindrical geometric space defect is applied to the electron be- 
haviour in the system which can be regarded as a simplified model of a double-wall nan- 
otube. By solving the Schrodinger equation in the region of space with cylindrical geomet- 
ric defect we explore the influence of such geometric defect on the energy gap and charge 
distribution. The effect is qualitatively similar to the one obtained earlier by means of 
traditional simulation methods. In general, the geometric approach can not compete with 
the known methods of theoretical study of the nanostructures, such as molecular dynamics. 
However it may be useful for better quahtative understanding of the electronic properties 
of the nanosystems. 

1 Introduction 

The theory of defects is an important part of the mathematical physics. In particular, 
topological defects attract much attention in modern condensed matter physics (see, e.g., 
[H 12] for an introduction and recent review). Another elegant approach in this area is the 
geometric theory of defects [3l H] (see also [5] for the introduction), which is formulated in 
terms of the notions originally developed in the theories of gravity. The main kinds of de- 
fects which were described in the framework of Riemann-Cartan geometry are dislocations 
and disclinations. This means that the curvature and torsion tensors are interpreted as 
surface densities of Frank and Burgers vectors and thus linked to the nonlinear, generally 
inelastic deformations of a solid. Recently, the qualitatively new kind of geometric defect 
corresponding to the cylindrical geometry has been described in [6]. In the framework of 



geometric approach, the static massive thin cyhndrical shells are considered as sources in 
Einstein equations. It turns out that these defects correspond to 6- and ^'-function type 
energy-momentum tensors. It would be interesting to find some useful applications of this 
new kind of defect. One may think about cosmological applications and/or about the 
condensed matter physics applications, which can be not infinitely far from each other [7]. 

One of the natural question to address is in which way the new kind of defects can be 
applied to the condensed matter physics? In this respect we can mention, for example, 
the recent works [SI IH] devoted to the cylindrical nanoparticles in a liquid crystal solvent. 
It would be interesting to see whether the cylindric defects can be useful in description of 
these systems. However, at the first place, the cylindrical geometry is associated with the 
nanotubes (see, e.g., [HI HOl [121 US] for the general review). As we shall see in what follows, 
the most appropriate object of our interest would be not the usual single-wall nanotubes but 
the double- wall nanotubes (DWNTs) or, more general, multi-wall nanotubes (MWNTs). 

The MWNTs attract great attention [HJ [15], [16] , in particular due to higher temperature 
stability and stiffness compared to the similar single-wall nanotubes. The inner and the 
outer layers of the DWNTs can be metallic (M) or semiconductor (S). Correspondingly, 
there are DWNTs of different kinds, namely: semiconductor-semiconductor (S/S case), 
metal-semiconductor (M/S case), semiconductor-metal (S/M case) and metal-metal (M/M 
case). It has been emphasized recently in [17] that the M/M type of the DWNTs is the 
most difficult to meet, but, at the same time, this particular configuration is especially 
interesting to explore. Let us note that there are some theoretical and experimental data 
available on the energy gap of a MWNTs compared to the purely metallic tubes [18] (see 
also [15]). 

In the present work we shall consider how the cylindrical geometric defects can be 
related to the theoretical investigation of the M/M type DWNTs. Our approach will 
be to consider the solution of the one-particle Schrodinger equation in the space with 
the cylindrical geometric defect. Let us note that the idea to explore the one-particle 
Schrodinger equation as a way to describe nanotubes is not completely new. One can 
mention, e.g., the work of Ref. [I9], where a nonrelativistic spinless particle moving in a 
cylindrical surface (also in a thick cylindrical hollow) is investigated. On the other side, 
there are also publications on quantized particles moving in a curved background, with 
or without defects. Many of them address the problem of a test particle in a spacetime 
with cosmic strings and are motivated by the corresponding cosmological models. Also, in 
Ref. [20] , for example, the Schrodinger picture description of vacuum states is studied and 
applied to simple cosmological models without cosmic strings. 

An important difference between the works mentioned above and our approach is that 
we do not consider the cylindrically symmetric potential. Instead, we consider the free 
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Schrodinger equation (without potential) in the space with the cyhndrical geometric defect. 
In order to see in which situation this method may be more natural, let us present the 
following observations. The conventional nanotubes are compounds where the electron 
can propagate only on the tube surface and not in the three-dimensional bulk. In the 
case of DWNTs there are two distinct conducting cylindrical shells. If we consider the 
electron in such system, the tube or space division shell between the two tubes is unlikely 
to create an electromagnetic potential difference between the two conducting layers while 
being essentiall for the electronic properties of the compound. Therefore, the presence of 
the division between the two shells fits the geometric defect case described in [6] , so it looks 
natural to perform some investigation of the corresponding quantum systems. 

It is necessary to mention some recent publications motivated by condensed matter ap- 
plications of geometric approach. In Ref. [2T], the effect of curvature is analyzed in particle 
scattering. Concerning the problems in the presence of defects, we can also mention, e.g., 
Ref. [22j, where a particle which moves in a magnetic field in a space with disclination 



and screw dislocation has been studied. Also, the paper [23j explores the bounded states 




of an electric dipole in the presence of a conical defect. Despite our work can be seen as 
continuation of this line, we consider a rather different approach to the geometrical defects 



Indeed, there are well-known methods for calculating electronic and mechanic properties 
for different compounds, 3D, 2D or ID nanosystems, in particular. For example, using 
the density functional (DFT) - based methods we can explore different systems (including 
periodical ones), and one obtains the desirable results with high precision [2l]. At the same 
time, it looks interesting to have an alternative, more simple (albeit potentially reliable) 
approach which could permit us to obtain important qualitative information, and maybe 
also help in qualitative understanding, e.g., of the electronic properties of the mentioned 
systems. In what follows we shall develop relatively simple and to great extent analytic 
method based on the geometric theory of defects. As we shall see in brief, this approach 
looks justified, for it enables one to arrive at the qualitative understanding of the origin 
of modifications of energy gaps and density distribution due to the presence of the defect 
shell between the conducting layers. 

The paper is organized as follows. In the next two sections we present a brief but to some 
extent pedagogical description of the tube dislocation in the linear elasticity theory and in 
the geometric theory of defects. Despite the content of these sections is essentially the same 
as the one of [6], we include it here com the sake of completeness. In section 4 we consider 
the Schrodinger equation and describe the results of its numerical solution. In particular, 
we compare the energy levels and charge density distributions of the MWNTs with the 
ones of the similar cylindrical tube without defect (it can be obtained analytically) and 
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Figure 1: Negative (a) and positive (6) tube dislocations. 

with the estabhshed earher features of the energy spectrum of the system [151 [18]. Finally, 
in section 5 we draw our conclusions. 

2 Tube dislocation in the linear elasticity theory 

Let us start by describing the tube dislocations in the simplest case of linear elasticity the- 
ory. Consider a homogeneous and isotropic elastic media as a three-dimensional Euclidean 
space with Cartesian coordinates x\ y^, where i = 1,2,3. The Euclidean metric is de- 
noted by 6ij = diag (+ -I- +). The basic variable in the elasticity theory is the displacement 
vector of a point in the elastic media, x G M.^. In the absence of external forces, 

Newton's and Hooke's laws reduce to three second order partial differential equations which 
describe the equilibrium state of elastic media (see, e.g., [25]). 

{l-2a)Aui + didjU^ = 0. (1) 

Here A is the Laplace operator and the dimensionless Poisson ratio cr (—1 < o" < 1/2) is 
defined as 

_ A 
2(A + /i)- 

A and n are called the Lame coefficients, they characterize the elastic properties of media. 

Raising and lowering of Latin indices i,j, . . . can be done by using the Euclidean metric, 
Sij, and its inverse, S''^ . Eq. ([T]) together with the corresponding boundary conditions 
enables one to establish the solution for the field u'^{x) in a unique way. 

Let us pose the problem for the tube dislocation shown in Figure [Ha. 



This tube dislocation can be produced as follows. We cut out the thick cylinder of media 
located between two parallel cylinders of radii ri and r2 (ri < with the 
the symmetry axis of both cylinders, move symmetrically both cutting surfaces one to the 
other and finally glue them. Due to circular and translational symmetries of the problem, 
in the equilibrium state the gluing surface is also the cylinder, of the radius which will 
to be found below. 

Within the procedure described above and shown in FiglUa we observe the negative 
tube dislocation because part of the media was removed. This corresponds to the case of 
ri < r2. However, the procedure can be applied in the opposite way by addition of extra 
media to as shown in FigJH^. In this case, we meet a positive tube dislocation and 
the inequality has an opposite sign, r\ > r2- It is important to note that the approach 
adopted here above concerns our treatment of the division of the shells of the DWNT as 
a geometric defect. Of course, it can not be seen as a practical prescription for building 
nanotubes of other objects. 

Let us start by calculating the radius of the equilibrium configuration r*. This problem 
is naturally formulated and solved in cylindrical coordinates r, ip, z. Let us denote the 
displacement field components in these coordinates by m'', m''', . In our case, m*^ = = 
due to the symmetry of the problem, so that the radial displacement field u'^{r) can be 
simply denoted as u'^{r) = u{r). 

The boundary conditions for the equilibrium tube dislocation are 



u\ „ = , u\ = 



du 



m 



r=0 ' lr=cx> ' 



dUp 



dr 



(2) 



The first two conditions are purely geometrical, and the third one means the equality of 
normal elastic forces inside and outside the gluing surface in the equilibrium state. The 
subscripts "in" and "ex" denote the displacement vector field inside and outside the gluing 
surface, respectively. 

Let us note that our definition of the displacement vector field follows [5], but differs 
slightly from the one used in many other references. In our notations, the point with 
coordinates y*, after elastic deformation, moves to the point with coordinates x*: 

y^^x'iy) = y' + u\x). (3) 

The displacement vector field is the difference between new and old coordinates, = 
X* — y*. Indeed, we are considering the components of the displacement vector field, ^^(x), 
as functions of the final state coordinates of media points, x*, while in other references 
they are functions of the initial coordinates, y*. The two approaches are equivalent in the 
absence of dislocations because both sets of coordinates x* and cover the entire Euclidean 
space M^. On the contrary, if dislocation is present, the final state coordinates x* cover the 
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whole while the initial state coordinates cover only part of the Euclidean space lying 
outside the thick cylinder which was removed. For this reason the final state coordinates 
represent the most useful choice here. 

The elasticity equations ([1]) can be easily solved for the case of tube dislocation under 
consideration. The Laplacian and the divergence in the cylindrical coordinates have the 
form 

1 1 12 

AUr = -dr{rdrUr) + ^d'^Ur + d'^Ur - - ^d^Uy, , (4) 

d,u' = ^driru'') + ^d^u'P + d,u' , (5) 

where the indices are lowered using the Euclidean metric in cylindrical coordinates, Ur = 
u^, u^p = u'^r'^ and Uz = u^- Let us note that the last two terms in 1^ are due to the 
geometric covariant nature of the Laplace operator, which is constucted on the basis of 
covariant derivatives. 

One can remember that only the radial component differs from zero. The angular ip and 
z components of equations ([1]) are identically satisfied, and the radial component reduces 
to the ordinary differential equation, 

-dJru)] =0, (6) 
.r J 



dr 

which has a general solution 



b 

u = ar 

r 



depending on the two arbitrary constants of integration a and b. Due to the first two 
boundary conditions ([2]), the solutions inside and outside the gluing surface are 

Uin = ar , a > 0, 

b ^ (7) 

Hex = , b> 0. 

r 

The signs of the integration constants correspond to the negative tube dislocation shown 
in FigHla. For positive tube dislocation, FigJH^, both integration constants have opposite 
signs, a < and 6 < 0. 

Using the solution (|2]) and the third boundary condition (|2]), one can determine the 
radius of the gluing surface, 

= -■ (8) 
a 

After simple algebra the integration constants can be expressed in terms of the initial 
radius of external and internal cylinders 

_ r^-ri _ J_ rl-rl _ In , . 
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where 

1 J r2 + ri 

I = r2-ri, and = — - — 

are the thickness of the removed cyhnder and the radius of the gluing surface, respectively. 
The first expression in (Q restricts the range of the integration constant, < |a| < 1. 

It is interesting to note that the above formulas are applicable for both negative and 
positive tube dislocations. One has to take / > and Z < in these cases, respectively. 
Indeed, the negative thickness in the last case means that, before the deformation, the 
external radius must be smaller than the internal one and that the proper deformation 
consists in inserting some matter between the cylindrical surfaces. In both cases we see 
that the gluing surface r = r^, lies exactly in the middle between the radii ri and r2. 

Finally, within the linear elasticity theory, Eq.(I7]) with the integration constants 
yields a complete solution for the tube dislocation. This solution is valid for small relative 
displacements, when l/ri <^ 1 and //r2 ^ 1. It is remarkable that the solution obtained 
in the framework of linear elasticity theory does not depend on the Poisson ratio of the 
media. In this sense, the tube dislocation is a purely geometric defect which does not feel 
the elastic properties. 

In order to use the geometric approach we have to present the results obtained above 
in other terms. For this end we compute the geometric quantities of the manifold corre- 
sponding to the tube dislocation. From the geometric point of view, the elastic deformation 
([3]) is a diffeomorphism between the given domains in the Euclidean space. The original 
elastic media M^, before the dislocation is made, is described by Cartesian coordinates 
with the Euclidean metric 6ij. An inverse diffeomorphism transformation x ^ y induces a 
nontrivial metric on M^, corresponding to the tube dislocation. In Cartesian coordinates 
this metric has the form 

dy^ dy^ 

We use curvilinear cylindrical coordinates for the tube dislocation and therefore it is useful 
to modify our notations. The indices in curvilinear coordinates in the Euclidean space 
will be denoted by Greek letters x^, fi = 1,2,3. Then the "induced" metric for the tube 
dislocation in cylindrical coordinates i£] 

where is the Euclidean metric written in cylindrical coordinates. We denote cylindrical 
coordinates of a point before the dislocation is made by where y without index 

stands for the radial coordinate and we take into account that the coordinates <y9 and z do 

iNotel 
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not change. Then the diffeomorphism is described by a single function relating old and 
new radial coordinates of a point y = r — u{r), where 

r , r < 

u{r)={ 5 (11) 

- - , r > . 

r 

It is easy to see that this function has a discontinuity ar* + 6/r* = / at the point of the cut. 
Therefore a special care must be taken in calculating the components of induced metric. 
It proves useful to introduce the function 

a , r < r^, , 

b (12) 

which is continuous on the cutting surface. This function differs from the derivative u'{r) 
of the displacement vector field u{r) defined in ( ITTj) by the 5-function 

u' = u'ij-) = v{r) — I5{r — r*) . (13) 

Now we define the induced metric for the tube dislocation as 

ds^ = (1 - vfdr^ + (r - ufd^^ + dz^ . (14) 

It is easy to check that this metric outside the cut agrees with the expression ( ITOj) . The 
difference between the two expressions is that ( JT4l) is defined on the cutting surface while 
(fTOl) is not. The volume element corresponding to (HM is 

\/\g\ = (1 — f)(r — u) , where g = detg^i, . 



Let us note that the metric (1141) differs from the one which results from the formal substi- 
tution of y = r — u{r) into the Euclidean metric ds"^ = dy"^ + y'^dip^ + dz^ by the square of 
the 5-function in the grr component. This procedure is required in the geometric theory of 
defects, because otherwise the Burgers vector can not be expressed as the surface integral 
[5]. Finally, the metric component grr{f) = (1 ~ ^)^ of tube dislocation is a continuous 
function, and the angular component (7^^ = (r — -u)^ has the jump at the surface of the 
cut. 

The displacement vector field (fT3ll and induced metric (IHj) represent the solution of 
the field equations of the linear elasticity theory ([T]) with the boundary conditions ([2]) 
describing tube dislocation. The solution does not depend on the elastic properties of the 
media due to the universality of the linear approximation. 
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3 Tube dislocation in the geometric theory of defects 



Our next step is to go beyond the linear approximation and consider the defects in elastic 
media with a tube dislocation. In the geometric theory of defects [U HI |271 [28] (see [5] 
for the for review) the defects in elastic media with a spin structure are the framework of 
differential geometry. The main assumption is that the elastic media is a three-dimensional 
manifold with a Riemann-Cartan geometry. This approach enables one to treat the non- 
linear regime and also enables one to deal with a continuous distribution of defects. The 
tube dislocation in the geometric theory of defects was introduced and considered in details 
in Ref. [6], so here we present only a short account of the formalism. 

The geometry is determined by the Einsteins equations with the source term given by 
the energy-momentum tensor. We remark that the formalism under discussion is diffeo- 
morphism invariant and thus, in contrast to the linear elasticity theory, the displacement 
vector field does not show up explicitly in the field equations. At the same time we require 
that the linear approximation gives standard result. Therefore a natural choice for the 
energy-momentum tensor is 



where the curvature terms are constructed by using the metric obtained within the elasticity 
theory, e.g., f|T4|) for the tube dislocation case. The straighforward calculations yeild the 
following relevant component of the above equation [6]: 



It is worth mentioning that the above calculation is non-trivial because the curvature terms 
contains ambiguous pieces, like 5^, but surprisingly they cancel in the final expression. 

It turns out that the solution in the geometric theory of defects in some appropriate 
gauge conditions coincides exactly with the metric flT^ . Thus, the result obtained within 
the hnear elasticity theory is quite general and we can adopt the corresponding geometry 
as a background for the quantum consideration. 

4 Schrodinger equation in the presence of tube defect 

Consider the quantum description of a spinless particle in the space with tube disloca- 
tion. The quantum properties of our interest are encoded into the time independent wave 
function, ip{r,(j),z) which is the solution of the stationary covariant Schrodinger equation. 




(15) 




(16) 



(17) 
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Here we use the notation x = y2rnE/h^ where m and E are mass and energy of the 
particle. A is the covariant Laplacian operator, given by 

M = -^d, iy^g^'d,^) . (18) 

Notice that the wave function behaves as a scalar field when subject to covariant derivative, 
because the particle is spinless. For this reason flTK]) boils down to eq. (jl]) in cylindrical 
coordinates and, furthermore, the torsion does not couple with the wave function [§. The 
straightforward calculations lead to the following form of Schrodinger equation in the space- 
time with the elastic deformations described in section 2, 

' ■ ^ ■ +Tr^ = -X'^. (19) 



(1 — i;)2 dr'^ (1 — vY \r — u 1 — v J dr (r — u)"^ dcj)'^ dz 

In order to obtain the energy spectrum and density distribution, we need to to solve the 
above equation separately in the two distinct regions, namely for r < r^, (this implies 
also y < ri) and for r > (in this case y > r2). We shall denote the corresponding 
solutions for the wave functions as ?/'™*(r) and ip'^^^lr). One has to take into account also 
the normalization and boundary conditions, along with the matching condition at r = r^., 
namely, 

V''"*(r,) = V^''"*(r,) , (20) 

[21) 



dr ) r=r-t \ dr 

It would be very difficult to solve the equation f|T9l) in a direct way, especially outside 
the tube defect. However the solution becomes rather straightforward if we perform the 
inverse coordinate transformation r y = r — u{r). In this case the Schrodinger equation 
becomes very simple and in fact it is just the one for the free particle, 

_|_ 1 f^V" _|_ 1 d'^'^P _|_ f^^"^ 2^ ^22) 
dy"^ y dy y'^ dcf)"^ dz"^ 

After the equation (!22|) is solved one has to perform the coordinate mapping, that means 
the direct elasticity transformation in order to get the solution of the equation ( !T9|) of our 
interest. It is important to keep in mind that the triviality of the equation f l22l) does not 
mean the triviality of the solutions of the equation of our interest (fT9l) . The whole point 
is that the solutions are coordinate dependent (while the equation is covariant) and the 
physical results correspond to the coordinates defined in eq. ([3]). 

2Note2 
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The solution of the free Schrodinger equation (!22l) can be easily obtained via the usual 
method of separating variables if we consider a finite cylinder between z = —L and z = L. 
We set ipdi, (p, z) = p{y)^{(f))Z{z) and arrive, after replacing this form into the equation 
and some standard algebra, at the solutions for and z dependences, 

$(0) = e'""^ and Z{z) = cos{kz) , (23) 

where m = 0,1,2,... and k = nn/L {n = 1,3,5,...). These solutions correspond to the 
boundary conditions Z{—L) = Z{L) = 0. The radial function satisfies the equation 

dy'^ ~^ y dy ) ^ ^ 

where = ~ k"^- The general solution for the last equation has the form 

p{y) = Mxy) + C2 YMy) , (25) 

where Jm and Ym are Bessel function of the first kind and of the second kind (or Neumann 
function) respectively, and Ci, C2 are arbitrary constants. The above solution can be 
mapped from the "old" coordinate y into physical coordinate r, as it was explained above, 
providing us an important information on the spectrum of the system and also on the 
charge density distribution. Let us discuss this in some details. 



4.1 Energy levels 

Let us consider the solution of Schrodinger equation for a particle confined in a thin region 
around the geometrical defect, say, in old coordinates, tq < y < ri, r2 < y < r^. The two 
separate solutions are 

p'^^'iy) = cMxy) + c2YUxy) , y<ri, (26) 
p'^^'iy) = c^J^ixy) + CiY^ixy) , y>r2- (27) 

After the necessary mapping, r = 7/ + we arrive at the exact solutions of the radial part 
of the equation flTI?]) . One should mention that a real finite-size tube does not have the 
symmetry along the z-axis, however one can consider a very long and thin tube and thus 
assume the translational z-symmetry as a good approximation. In our case this means an 
approximation of a large L, such that ~ 0. 

Consider the mapping, which was already discussed above, in the detailed form. In 
what follows, we shall drop the bars and write x instead of x- In order to write down the 
solution in terms of physical coordinate, we have to impose the boundary conditions in the 
form 

ClJ-niiX^o) + C2Ym{Xfo) = CsJmixn) + CiY^ixTs) = . (28) 
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These conditions give us, immediately, C2 = — ciJm(x?"o)/^m(x^o) and Ci = —c^Jmix'^^) I 
Ymix"^?^- The last relations must be inserted on the matching condition at r = r*. However, 
it turns out that the old matching conditions (12T]) are physically unacceptable, since they 
imply that the charge density, proportional to ^/g\1p\'^, has a discontinuity at the surface 
r^,. Thus, the following conditions would be more appropriate: 



(29) 



y=r2 



Here we have used the expression for the volume element in cylindric coordinates, \/ g{r) = 
r. 

Unfortunately, the above conditions provide only two relations between arbitrary (com- 
plex) constants, in contrast to the equations (12T]) . which provide four relations. Thus one 
should use another conditions (complex equations). The most natural choice is continuity 
of the covariant density ^/g\'^p\'^■ Another condition with direct physical interpretation is 
the current conservation, that is the continuity of the current y/gJ'^, that means J'^ = 0, 
where 

zn 



It turns out that the relations (l29l) or, in a more general form. 



g 



'l'g^^d,p{y) = g'/' g^" d,p{y) , (30) 



y=ri 



y=r2 



imply the above physical assumptions and provide four necessary relations. On the top of 
that we have the normalization condition 

= 1. (31) 

V 

The equations ( 130|) determine the discrete energy spectrum. Of course the results are 
strongly dependent on the values of parameters tq, ri, r2 and rs. In order to compare the 
results with the available data, let us us specify these parameters as follows: 

To = 10 nm, ri = 54 nm, r2 = 56nm and r3 = lOOnm. (32) 

Typically, the outer diameter of CNT ranges between 2 and 20 nm and inner diameter 
ranges between 1 and 3 nm, interlayer distance is 3.4 nm as it is established by the high- 
resolution TEM techniques. The values ( l32l) correspond to the typical spacing between 
the layers for nanotubes, confirmed, in particular, by the experiments using diffraction 
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Figure 2: An illustration of the general geometrical form of a MWNT with the inner 
diameter of 20 nm and outer diameter of 200 nm. 



and HRTEM image techniques [301 EI]- We consider here the possibility of describing 
MWNT, with several layers of carbon or other materials. Indeed, since we are looking 
for a qualitative correspondence with the results obtained by other methods, there is no 
much sense in performing calculations for large amount of possible nanotube diameters. 
However, since our method is really simple, such calculation can be easily performed, e.g., 
for any particular MWNT where the fast theoretical evaluation of the effect is needed. The 
system, as showed in Figure 121 can be treated geometrically, by identifying nanotubes with 
geometrical properties. In our consideration, a different kind of nanotube at r = 55 nm 
means a geometrical defect at the corresponding region. 

Let us note that, taking the radius magnitude close to the ones which are typically found 
in nanosystems, including MWNT's, we also assume that the (e.g. carbon) nanotubes layers 
are of the same chirality. Indeed, one can suppose that a nanotube located at the middle, 
at r = 55 nm, has different chirality, as well as, pehaps, few nanotubes in its vicinity [T6] . 



Let us find the energy spectrum for the above choice, which depends on the quantum 
number m. We shall start from the fundamental state, m = 0. One can write equation ( 130|) 
in terms of just one arbitrary constant, say, Ci, which can be fixed by using the condition 
(13T]) ). The corresponding equation is transcendental and rather cumbersome, so it is not 
convenient to write it here. The solution was performed with the help of the software 
Mathematica [32] . We have found the values for x = V'2mE / h which are solutions of eq. 
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(!30|) . The first three eigenvalues are placed in the left column of Table 1 . For the purpose 
of comparison, we put the corresponding eigenvalues for the simplest case of the condicting 
cylinder without geometrical defect, in the right column. Similar results are presented in 
Table 2 for the m = 1 case. 



X {nm ^) 


with defect {1 = 2 nm) 


without defect (/ = nm) 


Xi 


0.03383 


0.03314 


X2 


0.07016 


0.06858 


X3 


0.10612 


0.10377 



Table 1: Energy spectrum for m = with tube defect (left) and without defect (right). 
The quantities Xi) X2 and Xs correspond to the first three energy eigenvalues. 

From these results, one can see that the tube defect does shift, at the first place, the 
energy spectrum as a whole while it almost does not change the energy gap. Still there 
is a very small effect and the tendency is to reduce the gap slightly in all cases. Similar 
situation occurs for m = 1 (see Table 2). This shift depends on the size / = r2 — ri of the 
cutting region. In Table 3, we put the first three energy eigenvalues for the cases / = Onm, 
I = —2nm, I = 2nm and / = Anm, for m = 0. Of course, / = corresponds to the case 
without defect. The case I = —2nm corresponds to a negative tube dislocation (see Figure 
lb). 



X {nm ^) 


with defect {1 = 2 nm) 


without defect (/ = nm) 


Xi 


0.04009 


0.03941 


X2 


0.07490 


0.07331 


X3 


0.10982 


0.10748 



Table 2: Energy spectrum for m = 1 with tube defect (left) and without defect (right). 

It is remarkable that the energy levels do shift under the mapping from the "old" 
coordinates to the physical ones while one could naively expect that these levels remain 
the same. The reason is that (as we have already mentioned) the mapping does not change 
the general solutions of the Schrodinger equations, while our interest is concentrated on 
the particular solutions with the given boundary conditions. These conditions do change 
under the mapping and this makes the solutions of the equation (fT9ll really nontrivial. 
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Figure 3: Representation of the data of Table 3 in the plane x ^ ^- The first three eigenval- 
ues, Xi; X2 and xs are plotted for the cases / = —2 ram, / = Onm, / = 2nm and I = 4 ram. 
The dashed hues are supposed to describe the assumed linear dependence. 

Table 3: The three first energy levels for m = for the cases / = —2 nm, I = nm, 

I = 2 nm and / = 4 nm. 

How is the dependence of the energy levels on the size /? As suggested by the data of 
Tables 1, 2 and 3, positive / induces a slight increase of energy eigenvalues, while a negative 
/ lowers them. The exact function — 1)2,3) is unknown and all we can do is to 

plot the data from Table 3 as shown in Figure [3l This plotting clearly suggests a linear 
dependence between the energy levels x and the size /. 



One can notice, by straighforward computation, that the energy gaps, Axi := X2 — Xi 
and Ax2 '■= Xs ~ X2, are slightly increased for / > Onm and decreased for / < Onm. 
For example, for / = 4 nm, Ax2 is 4.4% greater than the corresponding value of Ax2 for 
/ = nm. The values of Axi,2 are shown in Table 4, and are plotted in Figure |H We can 
observe again a linear relation between Ax and /. Let us notice that, by assuming linear 
dependence, the straight line describing Axi(0 has a lower angular coefficient than the 
one corresponding to Ax2(0- 
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Figure 4: Representation of the data of Table 4 in the plane Ax x I. The two first energy 
gaps, Axi = X2 — Xi S'lid Ax2 = Xs ~ X2, are drawn for different values of /. The dashed 
lines are supposed to describe the assumed linear dependence. 



Table 4: The first two energy gaps, Axi = X2 — Xi and Ax2 = — X2, for m = for 

different values of I. 



4.2 Density Curves 

In the ideal gas approximation, the volume integral of the density distribution, 
/ y/9\'^{yj4'j ^)\'^ d(f)dzdr, is proportional to J ?/|p(?/)p (ir, with the constant of proportion- 
ality having dimensions of area. From the condition fj3T|) . |cip has dimensions of [length]"'^ 
and one can find (using Mathematica software) the following values for |ci|~^: 2265.29 
nm^, 8985.42 nm^ and 8151.17 nm^ for Xi, X2 and Xsy respectively, for m = 0, with 
tube defect. The corresponding values without defect, for xi, X2 and Xs are 2315.70 nm^, 
8017.33 nm^ and 11471.10 nm^, respectively. 

With these values, one can plot the normalized density against the coordinate y. In 
Figure El for the lowest energy eigenvalue, Xi, and m = 0, the normalized density is plotted 
as well as the normalized density without defect for comparison. Notice that the curve in 
the presence of defect has a cut because the point ri is identified with r2, so the interval 
ri < y < r2 is formally absent in the space with tube defect. The defect raises density in 
the central region. It is interesting to observe that this is not the unique effect caused by 
the defect. With the help of the Mathematica software, one can calculate the value of y for 



16 



20 40 60 80 100 

y (nm) 

Figure 5: Normalized density for m = and x = Xi- The continuous line has a cut and 
describes density in the presence of tube defect. The dashed line describes density without 
defect. 

which the density reaches its maximum value. This value is given by 2/max ~ 53.80 for 
the density without defect and ^max ~ 52.95 nm for the density with m = and / = 2nm. 
Thus, the defect not only raises the density distribution, but it also drags its peak to the 
left. 

Let us remark that the curves describing the particle density have not sjTumetry around 
the straight line ?/ = 55 nm. The peaks of the curves in Figure [5] are localized in the left 
half of the region 10 nm < y < 100 nm, and not in the middle of it. Actually, this property 
is essentially related to the properties of the Bessel functions. Indeed, let us remmember 
that the system is not unidimensional, such that it does not possess symmetry around the 
middle point. 



For X = X2 and still m = 0, the corresponding density curves are plotted in Figure |6l 
In the presence of the tube defect, the density is lower in central region but greater in the 
marginal regions, y ~ 30 nm and y ^ 80 nm. The shifts in the density distribution are 
small and of the same order in both cases, xi ^.nd X2, but in the last case the difference is 
richer. In the central region, the defect raises the density for x = Xi a-^id lowers the density 
for X = X2- Of course, for a larger /, these effects may become much more sensitive. 
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Figure 6: Normalized density for m = and X = X2- The continuous line has a cut and 
describes density in the presence of tube defect. The dashed line describes density without 
defect. 

5 Conclusions 

We considered the solution of the Schrodinger equation in the restricted region of space with 
cylindrical geometric defect. The region is confined between the two cylindrical shells and 
the diameter of the defect is intermediate between the ones of these shells. Geometrically 
this configuration resembles the double-wall nanotube (DWNT) and, therefore, it can be 
regarded as a simplified model of such nanotube. Due to the fact that the electron can 
freely move within the space between the two border shells, one can associate the system 
with the M/M-type DWNT. The solution of the Schrodinger equation shows the influence 
of the geometric defect on the energy gap and charge distribution. In particular, one 
meets a modified energy spectrum and the distribution of charge density, in the ideal gas 
approximation. These effects are qualitatively similar to the ones previously reported for 
some MWNTs in [18], where the calculations were based on the tight-binding method. In 
this respect we can conclude that the approach based on geometric defects, which does not 
take into account anything but the general form of the DWNT, is able to provide some 
(very restricted, indeed) information about the electron behaviour. 

It is obvious that the method based on geometric defects can not compete with the 
standard approaches based on molecular dynamics. The reason is that the geometric 
defects method can not take into account full details of the structure of the compound 
and is, in some sense, too general. At the same time this method may become much 
more interesting if one develops it further and, in particular, learns how to deal with more 



18 



sophisticated versions of geometric defects. In particular, it looks possible to take into 
account the chirality of the nanotube and, also, include the external magnetic field. We 
expect to consider these issues elsewhere. 
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